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Abstract 

A system of four coupled oscillators that exhibits 
unusual synchronization phenomena has been ana¬ 
lyzed. Existence of a one-way heteroclinic network, 
called heteroclinic ratchet, gives rise to uni-directional 
(de)synchronization between certain groups of cells. 
Moreover, we show that locking in frequency differ¬ 
ences occur when a small white noise is added to the 
dynamics of oscillators. 
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1 Introduction 

Phase oscillators are used as approxima¬ 
tions for the phase dynamics of coupled limit 
cycle oscillators in the case of weak cou¬ 
pling |Kuramoto, 1984[ 


Pikovsky et al., 2001 


Hoppensteadt and Izhikevic h, 1997) . They ex¬ 

hibit synchronization and clustering phenomena 
IKuramoto, 1984| Sakaguchi and Kuramoto, 1986) , 
even if coupling function consist of the first har¬ 
monic only. If the second harmonic of the coupling 
function is considered, it is possible to observe 
switchings between different clusterings as a result 
of an asymptotically stable robust heteroclinic cycle 
flHansel et al., 1993] , It is known that heteroclinic 
cycles are not structurally stable but they may exist 
robustly for coupled systems. This is due to the 
existence of robust invariant subspaces for certain cou¬ 
pling structures that may support robust heteroclinic 
connections that are saddle-to-sink on the invariant 
subspaces and form a heteroclinic cycles [ Krupa, 1997| 
|Ashwin and Field, 1999| Aguiar et al., 201 1| ~ Exis¬ 
tence of robust heteroclinic cycles or more generally 
heteroclinic networks in a system of three and four 


globally coupled phase oscillators have been analyzed 
in |Ashwin et al,, 20081 and in [jAshwin et al., 2006j , 
where an extreme sensitivity phenomenon to detuning 
of natural frequencies has been observed. Namely, os¬ 
cillators loose synchrony even for very small detuning 
of natural frequencies. |Karabacak and Ash win, 20T0) 
have considered the third harmonic of the coupling 
function and observed one-way heteroclinic networks, 
which are called heteroclinic ratchets. A heteroclinic 
ratchet is a heteroclinic network that, for some axis, 
contains trajectories winding in one direction only. 
Heteroclinic ratchets give rise to extreme sensitivity 
to detuning of certain sign. Namely, synchronization 
of a pair of oscillators is possible only when the 
natural frequency of a certain oscillator is larger than 
the other. We call this phenomenon uni-directional 
(de)synchronization. 

In the sequel, we identify a heteroclinic ratchet 
for a system of four coupled oscillators in Sec¬ 
tion [2] Although the system is less complicated 
than the original ratcheting system considered in 
]Karabacak and Ashwin, 20T0] , it exhibits more com¬ 
plicated dynamics: uni-directional synchronization be¬ 
tween groups of oscillators, explained in Sections[3]and 

m 

2 A Model of Four Coupled Oscillators That Sup¬ 
ports Heteroclinic Networks 

Consider the following well-known model of coupled 
phase oscillators: 


9 , = 


K \ 

TT ^ C, i Q ( O'i 0 t ). 


N 


0 ) 


Here, 9i € T = [0, 27r) denotes the phase of oscillator 
i and co, is its natural frequency. The connection matrix 
































{cij } represents the coupling. c t3 = 1 if oscillator i re¬ 
ceives an input from oscillator j and c l3 = 0 otherwise. 

Since g(-) is a 27r-periodic function, it can be 
written as a sum of Fourier harmonics: g(x) = 
Y^kLi r k sin (kx + a*). Without loss of generality, we 
may set K = N and rq = —1 by a scaling of time. 
Let us choose the following coupling function with two 
harmonics only: 


g(x) = — sin(x + ai) + r 2 sin( 2 x) ( 2 ) 

The model (QJ is used as the approximate phase dy¬ 
namics of weakly coupled limit cycle oscillators, and 
the weak coupling gives rise to a T 1 phase-shift sym¬ 
metry in the phase model (0. Hence, the dynamics of 
0 is invariant under the phase shift 

(0ii 02: • ■ •) On) > (0i + e, 02 +£,•■•, On + e) 

for any e € T. Below, this symmetry is used to reduce 
the dynamics to an (N — l)-dimensional phase differ¬ 
ence system. 

Let us consider the coupled phase oscillator system 
0 with the coupling structure given in Figure Q] This 
gives rise to the following system: 



Figure 1. An asymmetric coupled cell system: this gives the cou¬ 
pled system of form (3- 



Figure 2. Two balanced colorings of the coupled cell system in 
Figure |T] These give rise to 2-dimensional invariant subspaces X\ 
and X 2 of the system (4). 


0i = wj. + g(0 ± - 0 2 ) + g(0i - 0 3 ) 

02 = W 2 + <?(02 — 01 ) + < 7(02 — 0 < l ) 

03 = W3 +0(03 - 02) +0(0) 

04 = W 4 + 0(04 — 0l) + 0(04 — 82 ). 

Defining phase difference variables as 0i := 0 1 — 02 , 
02 := 02—04 and 03 = 03 — 04 , we obtain the following 
dynamical system for phase differences: 

01 = W13 + 0(01 + 03 - 02) + 0(01) 

-0(03 - 02) -0(0) 

02 = ^24 + 0( —01 — 03 + 02) + 0(02) , 

“S(-03 -0l) -0(-02) 

03 = W 34 + 0(03 - 0 2 ) + 0(0) 

-0(-03 -0l) -0(-0 2 )- 

u>ij denotes the detuning between oscillator i and os¬ 
cillator j, namely = Wj — 

We first assume identical oscillators, that is 

wi = • ■ • = W 4 = w => Wy = 0 Vi, j. (5) 

Oscillators with different natural frequencies will be 
considered in Section^ 


2.1 Invariant Subspaces 

The assumption that the oscillators are identical 
makes it possible to use the balanced coloring method 
[[Stewart et al., 2003) to obtain invariant subspaces of 
the system (j3j. A coloring of cells in a coupled cell 
system is called balanced if cells with identical color 
receives the same number of inputs from cells of any 
given color. A balanced coloring gives rise to an in¬ 
variant subspace obtained by assuming that the cells 
of same color have identical states. The converse of 
this statement is also true. Namely, for a given cou¬ 
pling structure, a coloring is balanced if the corre¬ 
sponding subspace is invariant under any system hav¬ 
ing that coupling structure. Therefore, the invariant 
subspaces obtained by the balanced coloring method 
are robust under the perturbations that preserve the cou¬ 
pling structure. For an introduction to this theory, see 
[Golubitsky and Stewart, 2006]. 

Using the balanced coloring method the invariant sub¬ 
spaces of the coupled cell system given in Figure[T|can 
be found as in Table Q] Using the above-mentioned 
phase difference reduction, the corresponding invari¬ 
ant subspace in T 5 for the system <|4j are also listed 
in Table Q] Note that for the system Q. there are only 
two 2-dimensional invariant subspaces, namely X\ and 
A+. Balanced colorings for these invariant subspaces 
are given in Figure [2] The invariant subspaces X\ , 
X 2 and their intersection A 3 can support a robust hete¬ 
roclinic cycle (see [Ashwin et al., 20 ill for robustness 
criteria of heteroclinic cycles). 















Balanced 

colorings 

Invariant subspaces of the system © 
on T 4 and system Q and on T 3 

{1|2|3|4} 

X 0 = T 4 

Ao = T 3 (whole space) 

{1|2|34} 

Ai = {9 G T 4 | d 2 = 6 i} 

Ai = {<j> G T 3 \(j >2 = 0} (0i — 03 plane) 

{12|3|4} 

X 2 = {9 G T 4 | Ox = 63 } 

X 2 = {0 G T 3 0i = 0} (02 — 03 plane) 

{12|34} 

X 3 ={6 GT 4 1 6 »i = 63,62 = 0 4 } 

A 3 = {0 G T 3 0i = 02 = 0} (03 axis) 

{134|2} 

X 4 = {9 G T 4 | 6 »i = 9s = e 4 } 

X 4 = {0 G T 3 0i = 03 = 0} (02 -axis) 

{1234} 

As = {6» G T 4 | 6 i = 6 2 = d 3 = 9 4 } 

As = {( 0 , 0 , 0 )} (origin) 


Table 1. Balanced colorings of the coupled system given in Fig¬ 
ure \T} The corresponding invariant subspaces for the system 0 and 
the corresponding reduced invariant subspaces for the system {4} are 
given. 



Figure 3. Phase portraits of the system s on invariant subspaces 
X\ (a) and X2 (b) are illustrated for parameters given in (8). Red 
lines indicate robust heteroclinic trajectories. Thick red lines are 
the winding heteroclinic trajectories. The robust heteroclinic ratchet 
formed by these winding and non-winding heteroclinic trajectories 
and the saddles p and q is shown in (c). 


2.2 Existence and Stability of a Heteroclinic 
Ratchet 

A heteroclinic ratchet (first defined in 
(Karabacak and Ashwin, 2010| ) is a heteroclinic 
network that contains a heteroclinic cycle winding in 
some direction and does not contain another hetero¬ 
clinic cycle winding in the opposite direction. To be 
precise, a heteroclinic cycle C C 1’ v parametrized 
by x(s) (x: [ 0 , 1 ) —> T w ) is winding in some di¬ 
rection if there is a projection map P: —> R. 

such that the parametrization x(s) (x: [ 0 , 1 ) R N ) 
of the lifted heteroclinic cycle C C satisfies 
lim^i P(x(s)) — P(x(0)) = 2 kit for some positive 
integer k. A heteroclinic cycle winding in the opposite 
direction would satisfy the same condition for a 
negative integer k (see (Ashwin and Karabacak, 201 1| 
for general properties of heteroclinic ratchets). 

As discussed above, the system © may have a robust 
heteroclinic network on the invariant subspaces X\ and 
X2 . Such a heteroclinic network should be connecting 
saddles on X 3 = X \ n X 2 . Reducing the equations in 
© to X 3 and considering identical natural frequencies, 
we get 


<j>3 = g(<h)-g{-<h)- ( 6 ) 

This system is Za-cquivariant, and therefore it can ad¬ 
mit a codimension -1 pitchfork bifurcation of the zero 
solution under some nondegeneracy conditions. The 
saddles emanating from this bifurcation are on A 3 
and they are of the form p = ( 0 , 0 ,j» 3 ) and q = 
(0,0 — p 3 ). The value P 3 can be obtained by solving 


d(P 3 ) ~ g(~P 3 ) = 0 as 


p 3 = cos 


-1 


cosai 
2r 2 cos a 2 


(7) 


In order to show that there exist heteroclinic connec¬ 
tions between p and q we use the simulation program 
XPPAUT (Ermentrout, 2002| . We identify a hetero¬ 
clinic ratchet for the parameter values 


<*1 = -5, r 2 = 0.15 (8) 

(see Figure [3]). On A|, the heteroclinic ratchet con¬ 
tains a non-winding trajectory and a trajectory winding 
along +(j>i and —0 3 directions (see Figure[3ji). On X 2 , 
it contains a non-winding trajectory and a trajectory 
winding along +0 2 direction (see Figure [3jt). These 
four connections and the saddles p and q form a hetero¬ 
clinic ratchet (see Figure^). For the parameters given 
in ®, p can be found as (0,0,0.3315). Considering 
the Jakobien of © at p, we can find the eigenvalues of 
the saddle p as A ( j P ' 1 = —0.3112, A ^ = 0.2967 and 
A^ = —0.0636. Similarly, the eigenvalues of q = 
—p = (0,0,0.3315) can be found as A ^ 9 " 1 = 0.3130 
and A^ } = -0.3276 and \ { 3 q) = -0.0636. These 
eigenvalues of p and q correspond to the eigenvectors 
Vi = (1,0,0) e Xi, v 2 = (0,1,0) G X 2 and u 3 = 
(0, 0,1) G A 3 . A heteroclinic cycle is attracting if the 
saddle quantity, defined as the absolute value of the ra¬ 
tio between the product of the eigenvalues correspond¬ 
ing to the expanding connections and the product of the 
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Figure 4 . A solution of the system 0 converging to the hetero¬ 
clinic ratchet. Initial states are chosen as ( 2 , 1 , 0 . 5 ). The solution 
shows the peculiar property for heteroclinic cycles that the residence 
time near equilibria increases as t —Y OO, before (f)\ and (f >2 get 
locked at zero due to the precision errors. 


eigenvalues corresponding to the contracting connec¬ 
tions is smaller than 1 [Melbourne, 19891 - Hence, we 
can conclude that the heteroclinic ratchet for the system 
© with parameters given in ® is asymptotically sta¬ 


ble since the saddle quantity 


\ (p) \ (•?) 


x(9) 


lambda, 


(p) 


= 0.9532 


is less than one. 

A solution of © converging to the heteroclinic ratchet 
can be seen in Figure [4] The increase in the residence 
time near equilibria is typical for a solution converging 
to a heteroclinic network. Winding of <j >i and (f >3 occur 
at the same time, respectively in positive and negative 
directions, due to the winding heteroclinic trajectory on 
X\ (see Figure^). Winding of rb-i occur in the positive 
direction due to the winding heteroclinic trajectory on 
X 2 (see Figure^). Since at each turn the solution gets 
closer to the equilibria p and q, after some time, cj)\ and 
02 get locked at zero due to the precision errors. Note 
that, the invariant subspaces Xi and X- 2 serve as bar¬ 
riers, and therefore no solution can pass through them. 
For this reason the solution winds in <j> 1 and 03 direc¬ 
tions only one time. Since winding in —<j >3 direction 
occurs together with the winding in +(f> 1 direction, this 
also happens only one time. However, these barriers 
can be broken by noise and/or detuning of natural fre¬ 
quencies leading to the uni-directional synchronization 
phenomenon. 


3 Uni-directional Synchronization in the Model 

We say that oscillators i and j are (frequency) syn¬ 
chronized if the observed frequency differences 


\of-of \ 

flij = lim -(9) 

£—>•00 t 

is equal to zero. Here Of £ R is the lifted phase vari¬ 
able for 6i € T. It is know that coupled oscillators can 
get frequency synchronized when the distance between 


Figure 5 . Uni-directional synchronization with respect to 6^13. The 
main graph shows the frequency differences 13, and £^34 for 
0 with parameters given in 0 as a function of detuning 0^13 when 
Cc?24 = CJ34 = 0 . Oscillators are frequency synchronized when 
^13 ^ 0 and the synchronization fails for oscillator pairs ( 13 ) and 
( 24 ) whenever CJ13 > 0 . The insets show time evolution of the 
phase differences (pi for a positive value of CJ13. 



Figure 6. Uni-directional synchronization with respect to 6 ^ 24 - The 
frequency differences H13, 1^24 and H34 for (0 with parameters 
given in 0 are shown as a function of detuning UJ 24 when UJ 13 = 
CC34 = 0 . Oscillators are frequency synchronized when u >24 0 0 
and the synchronization fails for the oscillator pair ( 24 ) whenever 
W24 > 0 . 


their natural frequencies, namely |wy| := |oji — ujj\ is 
small enough. If frequency synchronization of oscil¬ 
lators i and j occurs only when a specific one of the 
oscillators has greater natural frequency, namely for 
a certain sign of u>ij, we say that synchronization is 
uni-directional. Uni-directional synchronization phe¬ 
nomenon has been shown to occur for oscillator pairs 
when an asymptotically stable heteroclinic ratchet ex¬ 
ists in the phase space [Karabacak and Ashwin, 20T0} 
jAshwin and Karabacak, 20~TT) . 

For the system <j4]), we investigate the effect of detun¬ 
ings W 13 , W 24 and W 34 on the synchronization of oscil¬ 
lators, respectively in Figure [3 [ 6 ] and [3 Due to the 
winding connections in the heteroclinic ratchet, uni¬ 
directional synchronization occurs for detunings W 13 
and W 24 . However, because of the connection winding 
both in +<f> 1 and —(j )3 directions, a positive detuning 


















































Figure 7 . Bi-directional synchronization (the usual case) with re¬ 
spect to CJ34. The frequency differences ^13, ^24 and ^34 for 
0 with parameters given in 0 are shown as a function of detuning 
CJ34 when (jJ 13 = CJ24 = 0 - Oscillators are frequency synchro¬ 
nized when |co?341 is small enough and the synchronization fails for 
the oscillator pair ( 34 ) for large values of |Cc-734 1 . 


Detuning 

Direction 

Natural 

Frequencies 

Winding 

Direction 

Observed 

Frequencies 

+^13 

UJ, UJ, UJ 

+01 , — 03 

n+,n+,n,o+ 

+U >24 

UJ, OJ+ , UJ, UJ 

+02 


+U >34 

UJ+ , UJ, LJ+ , UJ 

0 

Q, Q, q, n 


Table 2 . The effect of detunings on the synchronization of oscil¬ 
lators for the system { 4 ). These can be obtained from Figures [ 5 ] [6] 
an dd Negative detunings have not been considered as they have no 
effect. UJ ^ ({l ) represents a number slightly larger than UJ ( 17 ). 

W 13 leads to synchronization of oscillators 1, 2 and 4. 
This is because 9\ — 0% = — ($3 — 64 ) ==> 9± = 6 * 4 . 
The synchronized groups of oscillators for each detun¬ 
ing case are given in Table 0 It is interesting that the 
oscillators 1,2 and 4 get synchronized for a positive 
detuning of 0143 , although the space {9 £ T 4 * * * * * 10 | 9\ = 
O 2 = 0 ,\} is not one of the synchronization spaces ob¬ 
tained by the balanced coloring method in Section I 2 TTI 
Hence, it is not an invariant subspace. 

4 Locking in Frequency Differences 

Noise induced uni-directional desynchro¬ 

nization of oscillators has been observed in 
flKarabacak and Ashwin, 20T0j for a coupled sys¬ 

tem admitting a heteroclinic ratchet. Here, we show 

that existence of a heteroclinic ratchet for the system 
i(4} leads to a locking in frequency differences when 

a small noise is applied. Figure [ 8 ] shows a solution 
of the system i[TJ under white noise with amplitude 

10 -12 . The noisy solution exhibits approximately 
equal number of windings in </>\ and &2 directions. 
This is because the noise is homogeneous and the 
invariant subspace X\ (resp. X 2 ) divides any e-ball 
around the equilibrium q (resp. p ) into two regions 
of attractions of equal volume for the winding and 



Figure 8. A solution of the system 0 for parameters given in 
0 with no detuning and with additive white noise (amplitude^ 
10 -12 ). 

non-winding trajectories. On the other hand, the num¬ 
ber of windings in 0 \ and --©3 directions are exactly 
the same because of the structure of the heteroclinic 
ratchet in Figure [3] 

As a result, the solution gives rise to the following 
frequency locking between frequency differences: 

fil 3 = SI 24 = ^34. (10) 

This is in agreement with the simulation results given 
in Figure[ 8 ] Therefore, the observed oscillator frequen¬ 
cies f \ := lim t \6 f(t)\/t are in the following form: 


(<h \ 


(n + s \ 

H2 


n 

H3 


n + s 

V^4 ) 


\n + 2Sj 


where 8 is a positive number. This type of a result can¬ 
not be seen directly from the connection structure of 
the coupled system, and is a consequence of the par¬ 
ticular heteroclinic ratchet that the system admits. Al¬ 
though the noise induces synchronization of oscillators 
1 and 3, the balanced coloring method explained in 
Section [2+1 does not give {9 € T 4 | 9\ = $ 3 } as an 
invariant subspace. 

5 Conclusion 

We have analyzed a system of four coupled phase os¬ 
cillators. The existence of an asymptotically stable het¬ 
eroclinic ratchet gives rise to uni-directional synchro¬ 
nization of certain groups of oscillators and induce a 
particular locking in the frequency differences of oscil¬ 
lators when small amplitude white noise is introduced 
to the system. These phenomena also lead to frequency 
synchronization of some oscillators, that can not be 
found by using the balanced coloring method, therefore 
does not correspond to any synchrony subspace. 

For future works, the relation between the connec¬ 
tion structure and possible synchronization groups can 











































be studied. Although the synchronization groups can 
not be inferred from the coupling structure directly, the 
coupling structure serves to create invariant subspaces 
on which heteroclinic ratchets can be supported. For 
this reason, the coupling structure plays an indirect role 
on the existence of possible synchronization groups. 
Another direction could be to study bifurcations of het¬ 
eroclinic ratchets that result in winding periodic orbits 
on torus. This can explain the effect of small detunings 
of natural frequencies on the observed frequencies in a 
complete way. 
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